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ABSTRACT 

We have analysed the magnitudes, kinematics and positions of a complete sample of 320 PNs in 
the elliptical galaxy NGC 4697. We show (i) that the PNs in NGC 4697 do not constitute a single 
population that is a fair tracer of the distribution of all stars. The radial velocity distributions, mean 
velocities, and dispersions of bright and faint subsamples differ with high statistical confidence, (ii) 
Using the combined data for PNs brighter than 26.2, we have identified a subpopulation of PNs which 
is azimuthally unmixed and kinematically peculiar, and which thus neither traces the distribution of 
all stars nor can it be in dynamical equilibrium with the galaxy potential, (iii) The planetary nebula 
luminosity functions (PNLF) of two kinematic subsamples in NGC 4697 differ with 99.7% confidence, 
ruling out a universal PNLF. We estimate that the inferred secondary PN population introduces an 
uncertainty in the bright cutoff magnitude of ~ 0.15 mag for this galaxy. - We argue that this 
secondary PN distribution may be associated with a younger, > 1 Gyr old stellar population, perhaps 
formed in tidal structures that have now fallen back onto the galaxy, as has previously been suggested 
for the X-ray point sources in this galaxy, or coming from a more recent merger/accretion with a red 
galaxy. The use of PNs for extragalactic distance determinations is not necessarily compromised, but 
their use as dynamical tracers of dark halos will require deep observations and careful analysis of large 
PN samples. 

Subject headings: elliptical galaxies:kincmatics- individual (NGC 4697)-galaxies: kinematics and 
dynamics- galaxies: distances and redshifts-planetary nebulae: general 



1. INTRODUCTION 

Planetary Nebulae (PNs) have become increasingly im- 
portant in extragalactic astronomy, for distanc e deter- 
mina ti ons via their luminosit y function (LF ) ( Jacobv 
| 1989f ICiardullo et all 119891 iJacobv Ciardullo Fordl 
ll990llMendezll999HFerrarese et alJ200(llCiardullol200.1 
and references ther ein), as kinematic trace r s of the dark 
halos of galaxies iArn aboldi et alJ Il99cl iSae lia et alJ 
1200(1 iMendez et all 1200 H iRomanowskv et all 1200^ 7 
and as tracers for the distribution and kinematics 
of the diffuse stellar population in galaxy clusters 
llFeldmeier. Ciardullo. fe Jacobv 1 119981 iFeldmeier et aD 
l2004HArnaboldi et al.M2002L 12004 IQerharH et 3,1.1 120051 . 
Due to their strong narrow line emission at [OIII]A5007, 
PNs can be easily detected out to distances be- 
yond 20 Mpc wit h narrow-band photometry and slit- 
jess spectroscopy llFeldmeier. Ciardu llo. fc Jacobv1ll998t 
Arnaboldi et alJl2002HMendez et alJl200lHDouglas et all 
20021) , an d to ~ 100 Mpc with multi-slit imaging spec- 
troscopy l)Oerhard et al.ll2005l) . Moreover, they are ob- 
served in elliptical and spiral galaxies, making them an 
indispensible tool to support distances obtained by other 
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methods (such as Cepheids, Surface brightness fluctua- 
tions, the Tully-Fishcr relation, SNe la), and to mea- 
sure the kinematics of stellar populations whose surface 
brightness is too faint for absorption line spectroscopy. 

For distance determination the Planetary Nebulae Lu- 
minosity Function (PNLF) is normally modeled as hav- 
ing a universal shape that depends only on the absolute 
bright magnitude cutoff M*: 

N(M) oc e °- 307M (l - e ^ M *-W) , (1) 

where N( M) is the number of P Ns with absolute mag- 
nitude M llCiardullo et alJll989|) . Observationally, the 
cutoff magnitude M* has a quasi-universal value of —4.5 
with only a weak dependence on host galaxy metallicity 
expressed by the system's oxygen abundance, and which 
can be compensated for by a quadratic re lation in [O/H] 
Ipopita et aU 119921 ICiardullo et alJl2002|) . In practice, 
the PN magnitudes m(5007), after correcting for the in- 
terstellar reddening, are fitted to the model PNLF of 
cq. ^ convolved with the photometric err or profile, yield- 
ing a value of the distance modulus ijCiardullo et alJ 
1989). The absence of any systematic variations in M* 
and the PNLF shape has been verified in galaxies with 
significant population gradients, and among galaxies of 
different morphologies within galaxy clu sters /groups up 
to Virgo l)Jacobvlll997l iCiardullol 120031 and references 
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therein) . 

This universality of the PNLF and the cutoff magni- 
tude M* must be considered surprising, given that the 
PN luminosity in the [OIII]A5007 line depends on the 
mass and metallicity of the central star, as well as on 
the electron gas temperature, optical thickness and dust 
extinction of the surrounding nebula. Indeed, some cur- 
rent semi-analytic simulations of the PN LF seem to be 
at od d s with the observational trends. IMendez et al.l 
1)1993(1 : IMendez fc Soffriel <|1997|) indicate small possible 
dependencies of M* on the total size of the PN popula- 
tion, on the time elapsed since the last episode of star 
formation, and on how optically thin the PNs are; con- 
cluding, however, that only careful studies would detect 
such effects in the observe d PNLF. In contras t, more 
recent PNLF simulations bv lMarigo et alJ 1)2004)1 contra- 
dict the observed narrow spread in M* and predict large 
variations of several magnitudes depending on a variety 
of realistic star formation and evolution scenarios. So is 
the PNLF truly quasi-universal and its cutoff magnitude 
nearly independent of population age and metallicity? 

PNs are also important as test particles to study the 
kinematics and dark matter distribution in the halos of 
elliptical galaxies. Since the PN population is expected 
to arise from the underlying galactic stellar distribu- 
tion, their radial velocities can be used as effective kine- 
matic tracers of the mass distribution. However, the re- 
quire d PN sample sizes are many 100s l|Merritt fc Sahal 
1993), or at least 100 or more in conjunction with ab- 
sorption line spectroscopy, which has li mited this appli- 
cation to only a few nearby galaxies <lHui et a ll |1995t 
Arnaboldi et alJll998l iSaglia et all 120001: IMendez et al l 
2001) IRomanowskv et all 120031 IPeng. Ford fc Freeman! 
2004). In recent simulations of disk galaxy mergers in- 
volving dark matter, stars, and gas. iDekel et alJ (2005) 
predict that the young stars formed in the merger have 
steeper density profiles and larger radial anisotropy than 
the old stars from the progenitor galaxies, and they ar- 
gue that if the PNs observed in elliptical galaxies were 
to correspond to the young population rather than to 
all stars in the simulations, their velocity dispersion 
profile would match the m easured dispersion profiles of 
IRomanowskv et all ()2003)1 . So do PNs really trace the 
stars and their kinematics in elliptical galaxies? 

Different stellar populations may have, and in general 
would have, different phase-space distributions in the 
same galaxy potential. The simplest approach for dy- 
namical modelling, taking the PN velocities as a random 
sampling of the stellar velocities, is however valid only 
when the PN population properties and their kinematics 
are uncorrelated. Except in special cases this also re- 
quires that the PNLF is independent of the stellar pop- 
ulation. Vice-versa, if there existed differences in the 
PNLF or the bright cutoff magnitude for different stellar 
populations, they would best be identified by studying 
the correlations between PN magnitudes and kinematics 
or positions of these tracers, in a single galaxy where all 
PNs are at the same distance. 

In this paper, we report on such a study in the ellipti- 
cal galaxy NGC 4697, an excellent target for this purpose 
because of the large sample of PN velocities known from 
IMendez et alJ l)2001f l. Our analysis shows the existence 
of distinct PN populations which differ in their kinemat- 
ics, brightnesses, and spatial distributions. This suggests 



that the answer to both the questions posed above may 
be 'no' - in general, different stellar populations may 
have slightly different PNLFs, and the observed PN pop- 
ulation in elliptical galaxies may not be a fair tracer of 
their stars. The paper is organised as follows: in §|3we 
review the properties and PN data of this galaxy and 
discuss the magnitude and velocity completeness of our 
sample. Our statistical analysis of these data is given in 
§ El where we demonstrate the inhomogeneity of the sam- 
ple in the space of velocities, magnitudes, and positions. 
S§ 14151 conclude our work, giving also a brief discussion 
of its implications. 

2. NGC 4697: PROPERTIES AND PN SAMPLE 

NGC 4697 is a normal, almost edge-on E4-5 galaxy 
located along the Virgo southern e xtension. From 
the m ulti-colour CCD photometry of iGoudfrooii et all 
(1994), the effective radius is R e fi = 95", the mean el- 
lipticity is 0.45, and the PA is constant, consistent with 
a near-axisymmetric luminosity distribution. Isophotal 
analysis shows that this galaxy has a positive a4 co- 
efficient su£gesting_j L _diskdike structure within 0.6R c ff 
(ICaTteHll9S7tlScorza fc BendeHlT995ll . 

From optical spectroscopy, its d ominant stellar pop - 
ulation has an age of - 9 Gyr l)Trager et all ETOOOL 
consistent with the red mean B-V=0.91 colour. Stel- 
lar absorption line kinematics along the major axis 
(PA 66°) of NGC 4697 have bee n reported by 
iBinnev. Davies fc Illingworthl 1)1990)1 and lDeionghe et all 
()1996)1 : these velocity data can be well described by dy- 
namical models based on the luminous mass distribution 
only. 

IMendez et alJ 1)2001)1 detected and measured radial ve- 
locities of 531 PNs extending out to 300" (~ 3R e fi) in this 
galaxy, with observational errors of ~ 40kms _1 . Via 
dynamical analysis, they determined a constant mass to 
light ratio Tb = 11 within ~ 300", which is consistent 
with a 10 Gyr old stellar population with a Salpeter mass 
function and slightly super solar metallicity. 

X-ray observations with ROSAT 

(Sansom, Hibbard & Schweizer 2000) show only small 
amounts of hot gas in the halo of this galaxy. Using 
more recent Chandra data, iSarazin. Irwin fc Bregmanl 
(2000) could resolve most of the X-ray emission into 
nommiformly distributed X-ray binary point sources 
(XPS), suggesting that NGC 4697 has lost most of 
its interstellar gas. Though NGC 4 697 does no t show 
any signature of recent interactions, iZezas et"aTI ( 2003) 
present evidence that the distribution of these X-ray 
sources is inconsistent with the optical morphology of 
NGC 4697, and propose that these sources are mostly 
high mass X-ray binaries (HMXBs) associated with 
young stellar populations related to fallback of material 
in tidal tails onto a relaxed merger remnant, or to 
shock-induced star formation along tidal tails. They 
estimate typical fallback times of such tidal tails to be 
much longer than the settling timescale of the remnant 
and expect similar results for other elliptical galaxies 
with populations of ~ 10 Gyr age. 

For the work in this paper w e use the PN sample pre- 
sented by Mendez et al. (2001). After removing the pos- 
sible contaminants and unclear detections, they report 
unambigous detection of 535 PNs. However only 526 
out of 535 PNs have confirmed measurements of velocity 
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and magnitude, and we use these in our analysis. In or- 
der to determine the PNLF, it is crucial to estimate the 
magnitude where PN detection incompleteness sets in. 
Detectability of a PN varies with the background galaxy 
surface brightness; for a statistically complete sample the 
surface numbe r density of PNs sho uld be directly pro- 
portional to it. IMendez et all (|2001f) show that their PN 
sample is statistically complete down to m(5007) = 27.6 
magnitudes outside an elliptical region of semi-major axis 
60". In our analysis, we have thus defined two data sets: 
a Complete sample (with PNs brighter than 27.6, outside 
the central ellipse of semi-major axis 60"), and a Total 
sample (consisting of all PNs with measured magnitude 
and radial velocity). The total number of PNs in these 
data sets is 320 and 526, respectively. 

The systemic velocity (V sys ) of NGC 4697, obtained 
by averaging the observed velocity of all 526 PNs is 
1274 km s - 1 , which agrees with the values quoted in the 
literature iMendez et al. 1 120011 and references therein). 
The on-band filter configuration used to detect and mea- 
sure velocities of these PNs has a peak wavelength of 
5028 A, peak transmission of 0.76, equivalent wi dth of 
48.5 A, and FWHM of 60 A ijMendez et alJ l2001). The 
FWHM corresponds to a velocity range of ±1800kms _1 
around the systemic velocity of NGC 4697. Thus the 
filter transmission width is large enough to facilitate ob- 
servations of PNs with all velocities bound to NGC 4697, 
irrespective of their magnitude. Indeed, even at magni- 
tudes as faint as m(5007) = 28 in the Total sample, PNs 
with large velocities ~ 300 km s _1 are detected. Thus the 
velocity coverage in both samples (Total and Complete) 
is not biased with respect to the PN mag nitudes. 

The PN magnitudes were measured bv IMendez et alJ 
(2001) from their undispersed images; they are accurate 
to 0.1 and 0.2 mag for m(5007) brighter and fainter than 
26.5, with systematic effects b elow 2%. As a furthe r test 
relevant for the present work, IMendez et alJ l)2001l) used 
the redundancy provided between their E and W fields: 
plotting magnitude differences between the two measure- 
ments (E and W) of PN candidates as a function of dif- 
ference in distance from the center of the CCD, they 
found a scatter diagram wit hout any evidence of corre- 
lation. IMendez et alJ l)2001l) estimated the errors in the 
PN velocities from calibration, image registration, spec- 
trograph deformation and guiding errors to be 40 km s _1 . 
The velocities of 165 PNs were measured i ndependently 
in the E and W field of lMendez et alJl|2001fi : these veloci- 
ties agree within a standard deviation of 36 km s _1 . In or- 
der to check whether a systematic difference between the 
velocities of bright and faint PNs could be introduced by 
an asymmetric PSF (a possibility suggested by K. Free- 
man), we have superposed the PSF's of three groups of 
10 of the brightest PNs, one selected at random, and 
two selected among those PNs with the highest and low- 
est radial velocities. In the three cases we estimated the 
shift of the centroid of the entire PSF with respect to the 
centroid of the upper part. The shifts were smaller than 
lOkrns -1 , and in some cases they were in the opposite 
sense compared to the results discussed below. 

3. THE DISTRIBUTION OF NGC 4697 PNS IN VELOCITY, 
MAGNITUDE, AND POSITION 

In this Section, we search for stellar population effects 
in the kinematics of the PNs in NGC 4697, by analysing 



the Total and Complete data sets with respect to their 
three observables: velocity, magnitude and position. For 
both data sets, we convert the observed PN radial ve- 
locities into co-rotating or counter-rotating velocities, as 
follows. With the galaxy center at the origin of the ref- 
erence frame, and the X-axis oriented along the major 
axis (PA= 66 deg), the absorption line stellar- kinematic 
data predict positive line-of-sight mean velocity with re- 
spect to the galaxy systemic velocity, at slit positions 
towards the South- West of the center with X coordinate 
> 0, and vice versa. We denote this sense of rotation as 
co-rotating, and the opposite sense as counter-rotating. 
By definition, the major axis absorption line data is co- 
rotating. 

3.1. Sub-populations of PNs in the velocity-magnitude 

plane 

After subtracting the systemic velocity from the 
PN radial velocities, we define reduced velocities 
U = (V — V sys ) * sign(x) and denote the PNs with U > 
(<) as co- rotating (counter-rotating). By definition, 
the major axis absorption line data have U > 0. The 
resulting values of U are displayed against the observed 
magnitudes in Figure^ Even at magnitudes as faint as 
m(5007) = 28 in the Total sample, PNs with large ve- 
locities ~ 300 km s" 1 are detected, showing that there is 
no kinematic bias at faint magnitudes. Henceforth, un- 
less stated otherwise, we will always use the Complete 
sample for our analysis. 

Figure shows that the Complete PN sample appears 
to exhibit a correlation between magnitudes and kine- 
matics, with faint PNs showing more co-rotation than 
bright PNs. We have performed several statistical tests 
to verify the significance and look for the origin of this 
correlation. Table ^ shows the results of Pearson's r- 
test for correlated data. Velocities of counter-rotating 
PNs are strongly linearly correlated with their bright- 
ness, while those of co-rotating PNs are almost indepen- 
dent of their magnitude distribution. Further, we di- 
vided our sample into 3 equal magnitude bins each of 
size Am = 0.7, hereafter referred to as faintest, inter- 
mediate, and brightest PNs, and computed the mean re- 
duced velocity and its dispersion in each of these bins 
along with the significance of their differences. As shown 
in Table 13 the faint and bright PN subsamples defined 
through these bins have different mean reduced velocity 
and dispersions at 94% and 99.3% confidence, respec- 
tively. In Figure |21 we have plotted the cumulative ve- 
locity distribution of PNs in the brightest and faintest 
magnitude bins. There is a visible excess of bright PNs 
with counter rotating velocities. It is particularly evident 
from this figure that the brightest counter rotating PNs 
display a velocity distribution that differs from the rest of 
the PNs from the Complete sample with high confidence. 

Thus it is clear that the observed correlation between 
the PNs kinematics and their magnitudes is compelling, 
and it arises because the faintest and brightest PNs 
have significantly different velocity distributions. There 
appears to exist an additional component of bright, 
counter-rotating PNs with respect to the overall sample. 

3.2. Spatial Distribution 

If these correlations have a physical origin, they should 
also be manifest in the spatial distribution of these 
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Fig. 1. — PN reduced velocity U versus magnitude plots for the Total (T) and Complete (C) samples. Reduced velocities have been 
folded such that positive (negative) velocity denotes co- (counter-) rotating PNs. Red stars (with error bars) are mean U values of PNs 
binned in magnitude such that each bin contains 30 (50) data points in the Complete (Total) sample. The error bars include the sample 
variance in the bin and the velocity errors (40kms ) added in quadrature. The black solid lines are least square linear fits to this binned 
data. In both samples, PNs with m(5007) < 26.5 have mean U ~ 0. The dotted vertical lines separate the Complete sample into 3 bins, 
each of size 0.7 magnitude (see text). The Complete data set shows a linear correlation between m(5007) and U: Fainter PNs show more 
co-rotation than their brighter counterparts. In the Total sample, the vertical solid line at m(5007) = 27.6 denotes the limiting magnitude, 
though not all PNs brighter than this magnitude constitute the Complete sample. 

any kinematic evidence (like a rotation curve signature) 
relating these bright PNs to the central stellar disk: ei- 
ther they are not physically related to the disk, or the 
evidence from the data is inconclusive. 

Outside the incompleteness ellipse, the distribution of 
bright PNs does not follow the surface brightness of 
the host galaxy: there is a significant left-right asym- 
metry, with more bright PNs to the right side of the 
galaxy minor- a xis (2 7 ± 5) than to the left side (15 ± 4). 
iMendez et all (|200lD discuss at length the possible dif- 
ferences in their E (east) and W (west) fields, and are 
convinced that the maximum systematic errors in the 
measured photometry, positions and radial velocities are 
below 2%, 0".06 and 20kms~ 1 , respectively. Hence we 
conclude that the left-right asymmetry in number counts 
of bright PNs is not affected by detection uncertainties. 

Subsequently, we carried out several tests to check 
whether the brightest and faintest, or the co- and 
counter-rotating PNs are distributed differently in the 
galaxy. It turns out that the radial PN distribution is 
independent of their sense of rotation. However, the dis- 
tribution of PN distances from the galaxy mid-plane dif- 
fers for co- and counter-rotating PNs at a confidence level 
of 73%. 

The left-right asymmetry is confirmed by inspect- 
ing the azimuthal distribution of the faint and bright 
PNs. In the liter ature we found a related analysis by 
iZezas et alJ l|2003|) who compared the azimuthal distri- 
bution of Chandra X-ray point sources (XPS) with the 
optical surface brightness of NGC 4697. We follow their 
PA convention, and plot the cumulative angular distri- 
butions of the bright, intermediate, and faint PNs in our 
Complete sample in Figure [S] For compa rison, the right 
panel of Fig. 2 from IZezas et al.l ((2003) is also shown. 
The angular distribution of all PNs in the Complete sam- 
ple has a shape somewhere in between that of the XPS 
and that of the optical light. The brightest PNs are in 



TABLE 1 

Pearson's r-test for linear correlation of PN magnitudes 
with reduced velocity u, for the entire complete sample, 
the co-rotating, and the counter- rotating subsamples. 
Values of r close to ±1 indicate a strong linear 
correlation; values close to indicate little or no 
correlation. vr is the probability that two uncorrelated 

variables would give the r-coefficient as large as or 
larger than the measured one, for a normal distribution 
of r. Small values of P t imply significant correlation. 



PNs. Thus we now enquire whether the PN kinemat- 
ics and magnitudes depend on their spatial location in 
the galaxy. 

In Figure we plot the spatial locations of all the 526 
PNs in this galaxy. The central incompleteness ellipse 
is also displayed. PNs brighter than m = 26.2 (which is 
0.71 magnitude deeper than the brightest PN) are shown 
as filled blue squares (co-rotating) and filled red triangles 
(counter-rotating). Inside the incompleteness ellipse, the 
distribution of the bright PNs appears to be concentrated 
around an elliptical annulus. However, we did not find 
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TABLE 2 

Mean reduced velocity (U) and dispersion (a) of brightest 

AND FAINTEST PNS, ALONG WITH THE RESULTS FROM STUDENT 
T-TEST AND F-TEST. VERY SMALL VALUES OF Vt AND Vf IMPLY 
THAT THE DIFFERENCES BETWEEN THE OBSERVED U AND 0~ ARE 
STATISTICALLY SIGNIFICANT. 



PN sub-populations in NGC 4697 



5 



E <*• 



c 




3?n 

. 42 

169 


V 
f 


KS F & B: 0.02 


7 


/ 




' / 
1 Jf 

1 6 - 
'' f' 




> / f 
r S > 
/ J 









-400 -200 



200 400 




-100 -200 -300 -400 100 200 300 400 
U (km s~ 1 ) U (km s" 1 ) 




400 200 -200 -400 

U (km s"') 

Fig. 2. — Cumulative velocity distribution of brightest PNs (m < 
26.2, orange dashed line), faintest PNs (m > 26.9, purple dashed- 
dotted line), and the entire Complete sample (black solid line). The 
entire velocity range is shown in the top panel. In the middle panel 
the velocities are divided according to their sense of rotation. The 
bottom panel shows the cumulative distribution (shown from +ve 
to —ve values), normalised at U = 0. The Kolmogorov Smirnov 
probabilities show that the brightest counter rotating PNs have 
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ange line), and the faint subsample (purple dash-dotted line). The 
KS probability shows that these are all consistent with one another. 

complete disagreement with either of these distributions; 
they seem to be more concentrated in a narrow angular 
sector between 110° < PA < 150° (see Fig.©, with only 
1.9% probability that the faint and bright PN subsam- 
ples are drawn from the same azimuthal distribution. At 
the same time, the radial distribution of the faint and 
bright subsamples are not significantly different (Figure 

n. 

The velocity distribution of the brightest PNs is also 
correlated with their azimuthal distribution. In FigureEl 
we plot the mean radial velocity and its dispersion in an- 
gular sectors containing approximately constant number 
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Fig. 5. — Azimuthal distribution of discrete sources in NGC 
4697. Top: Cumulative distribution of azimut hs of X-ray point 
sources (black solid line), adapted from Zezas et al.l |20PS). The 
green dotted line is for a uniform distribution, while the red dashed 
line shows the azimuthal distribution of the galaxy optical light 
(in R-band). Bottom: PN s in the brightest, intermediate, and 
faintest bins as defined in § 13.11 are shown by the orange dashed 
line, the red dash-triple dotted line, and the purple dash-dotted 
line, respectively. All PNs in the Complete sample are shown by 
the solid black line. Note the large increase in the distribution 
of bright PNs and the moderate increase in the distribution of 
intermediate PNs, relative to the faint PNs, at PA ~ 110° (see text 
for details). The KS probability of the faintest and brightest PNs 
being drawn from the same underlying distribution is only 1.9%, 
while the intermediate PNs are still compatible with the faint PNs 
(KS probability 28%). 

of PNs from the bright subsample, as well as an angular 
running average of their mean radial velocity. Along both 
sides of the major-axis of the galaxy, the brightest PNs 
have positive velocity, while showing a relatively low ve- 
locity disper sion, perhap s due t o infall as suggested for 
the XPS bv iZezas et al.l (|2003). Along the minor-axis 
they have a large dispersion with a mean velocity that is 
negative. This kinematics is compatible with neither the 
faint PN velocities nor the stellar absorption line data. 
On the other hand, the fainter PNs show a regular az- 
imuthal distribution in the mean line-of-sight radial ve- 
locity and velocity dispersion. Note that their velocity 
dispersion is in the mean smaller than that of the bright 
sample. The kinematics of the intermediate luminosity 
PNs is intermediate between those of the bright and faint 
subsamples. The intermediate luminosity PNs thus con- 
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Fig. 6. — Azimuthal distributions of PN velocities. The PA is 
measu red clockwise from North, as in the convention of Zez as et alj 
(2003), so that the galaxy major axis is along PA = 114°, 294°. 
Solid (dashed) histograms show the mean velocity (dispersion) and 
their errors in different angular bins. The red ragged lines show 
running averages of mean velocity. Top: Velocities of PNs brighter 
than m(5007) = 26.2 in the Complete sample. Bin sizes are chosen 
so as to have approximately constant number of PNs (~ 7) in each 
bin. The mean velocity along the major (minor) axis is positive 
(negative) on both sides of the center. The running average is also 
over 7 PN velocities. The velocity dispersion is largest (smallest) 
on the minor (major) axis. Middle: Velocities of PNs with 26.2 < 
m(5007) < 26.9 in the same angular bins as in the top panel. The 
running average is over 27 PNs. Bottom: Velocities of PNs with 
26.9 < m(5007) < 27.6 in the same angular bins as in the top 
panel. The faint PNs show a rotation pattern as for the absorption 
line data but with a smaller peak velocity, and are consistent with 
a flat dispersion of ~ 120 kms -1 . The running average is over 33 
PNs. The blue dotted line shows a sinusoidal fit. The kinematics of 
the intermediate brightness PNs is intermediate between the faint 
and bright PNs. 
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Fig. 7. — Velocity within ±30" of the major axis for bright 
and faint/intermediate PNs separated at m(5007) = 26.2. Blue 
circles show stellar absorptio n line spectroscopy (ALS) data from 
Binncy, Davics & Illingworth 1 1990), red crosses are individual PN 
velocities, and grey diamonds are near neighbour running averages. 
The faint/intermediate sample is in approximate agreement with 
the ALS data, while the bright sample has a velocity consistent 
with zero within R c ff and positive velocities on both sides of the 
galaxy center. 

tain significant contributions from both of the different 
populations that dominate the bright and faint bins, re- 
spectively. 

Furthermore, including the bright PNs in deriving 
kinematics for the whole PN sample introduces signfi- 
cant contamination effects. In Figure we demonstrate 
this by plotting the mean PN velocities along a ±30" 
wide slit about the major-axis. The PNs have been sep- 
arated into faint /intermediate and bright sub-samples at 
m(5007) = 26.2. Inside the R e fT, the mean velocities 
of the faint PNs agree wit h stellar absorption line spec- 
trosco pic (ALS) data from lBinnev. Davies fc Illingworthl 
(1990), while the bright PNs have much smaller mean 
velocities, consistent with ~ 0. In the outer parts, the 
velocity asymmetry in the bright PNs leads to a posi- 
tive mean velocity on both sides of the galaxy center. 
Both the faint/intermediate sample and the entire PN 
sample also show some asymmetry: the outer mean ve- 
locities on the X < side reach zero, but not the neg- 
ative values expected from reflecting the positive values 
at X > 0. This confirms that also the intermediate/faint 
subsamplc contains a fraction of PNs stemming from the 
out-of-equilibrium population traced by the bright PNs. 



However, the streaming velocity of the majority of the 
faint PNs does appear to decrease on both sides of the 
center. Similarly contaminated results could be expected 
for the derived velocity dispersions. 

Several important conclusions can be drawn from these 
figures, (i) The bright PNs as defined in Section ITT7T1 and 
Fig- H do not trace the azimuthal distribution of light in 
NGC 4697. (ii) They do not trace the fainter PNs in their 
azimuthal kinematics; thus, a large fraction of them must 
belong to a separate PN subcomponent originating from 
a separate stellar population, (iii) Third, they are not 
in dynamical equilibrium in the gravitational potential 
of NGC 4697. (iv) Because this subpopulation does not 
trace the stars, including its PN velocities into dynamical 
analysis of the galaxy will lead to significant errors in the 
results. 



3.3. Kinematic dependence of PNLF 

What can we learn about the luminosity functions of 
the two kinematic components of the PN system in NGC 
4697? The cumulative magnitude distributions of the co- 
rotating and counter- rotating subsamples in Fig. ^ have 
only 28% KS probability of stemming from the same dis- 
tribution (top panel of Figure|SJ). However, the fact that 
there are both co-rotating and counter-rotating bright 
PNs in the bright sample whose azimuthal distribution 
is unmixed, shows that counter-rotation is not a clean 
discriminator after all for the secondary population of 
PNs in NGC 4697. Thus the luminosity functions of the 
main and secondary PN populations in this galaxy may 
be a lot more different than this figure of 28% would 
suggest. 

From Fig. we estimate that the main population of 
PNs in NGC 4697 has a mean radial velocity t7f a mt — 
45 kms -1 sin(PA — 24°), with a dispersion of erf a i nt ~ 
120 kms -1 , so its reduced mean velocity is < J7f a ; n t J$ 
45 kms -1 . Thus in the bottom panel of Fig. we show 
the cumulative magnitude distribution of the PNs in the 
velocity range — 100 kms -1 <U< 200 kms -1 and com- 
pare it with the cumulative distribution of the PNs with 
U < — 100 kms -1 . These velocity ranges are dominated 
by the main and secondary PN populations in the sam- 
ple, respectively. Now the Kolmogorov Smirnov signifi- 
cance test shows that the magnitude distributions from 
these sections of Fig. ^ have only probability 0.3% of 
being drawn from the same distribution. This result is 
strong enough to imply that the PNLF cannot be uni- 
versal - the PNLF in NGC 4697 depends on a kinematic 
selection. 

In the following, we will use the velocity range 
— 100 kms -1 < U < 200 kms -1 as an approximate kine- 
matic selection criterion for the main PN population 
in NGC 4697. The luminosity function of the strongly 
counter-rotating PNs in Fig. [I]- a first approximation to 
the luminosity function of the secondary PN population 
in NGC 4697 - differs from that of the main PN distribu- 
tion so defined in the sense that it contains more bright 
PNs near the cut-off and fewer faint PNs than the main 
population (see Fig. [SJ. Now an important question is: 
in what proportion do both populations contribute to the 
brightest PNs, and do they have different cutoff magni- 
tudes? 

To investigate this, we show in Figure [5] the veloci- 
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Fig. 8. — Top: Cumulative distributions of PNs magnitudes for 
all (black solid line), co-rotating (blue dotted line) and counter- 
rotating (red dash-dotted line) PNs in NGC 4697. The magnitude 
distribution of co- and counter-rotating PNs is similar with only 
28% probability. Bottom: Cumulative magnitude distributions of 
the main PN sample (-100 km s -1 < U < 200 km s -1 , blue dotted 
line) and the extreme counter-rotating sample (U < — lOOkrns - 1 , 
red dash-dotted line) defined in the text, along with that of the 
total sample (black solid line). The first two luminosity functions 
are different with 99.7% confidence; thus the PNLF cannot be uni- 
versal. 

ties, magnitudes, and position angles of the entire bright 
PN subsample (see Figs. ^ and ®. The top panel shows 
(i) that 13/16 of the bright PNs, whose radial velocities 
differ most from those of the faint population, are coun- 
terrotating. This explains the differences between the 
velocity distributions of co- and counter-rotating PNs 
that first suggested more than one PN population in 
Section O (ii) Also, even in the kinematically normal 
bright PNs, there is a large overdensity (10/12) in the an- 
gular range PA= 100° — 150°. The bottom panel shows 
in addition that many of the brightest PNs are either 
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Fig. 9. — Disentangling the bright PN subsample from ra- 
dial velocities, luminosities, and position angles on the sky. Top: 
Radial velocity - position angle plane. Red squares (blue dia- 
monds) show PNs whose reduced velocities are outside (within) 
— lOOkrns -1 < U < 200 kms -1 . Circled symbols denote counter- 
rotating PNs; many of these are kinematic outliers. The two grey 
lines denote the position angle of the minor axis. The distribution 
of blue diamonds is approximativcly uniform in PA except near 
PA= 100° — 150° where there are 12 PNs (crossed diamonds) in- 
stead of the expected 2 PNs. Bottom: Radial velocity — magnitude 
plane. The symbols are as in the top panel. 6/8 of the brightest 
PNs are either kinematic outliers or are found in the overdense 
angular region. 

kinematic outlyers or found in the angular overdensity 
(6/8). 

Clearly, to arrive at a main population of NGC 4697 
PNs that is in dynamical equlibrium in the gravitational 
potential, we must remove the angular overdensity. Then 
we are left with 16 PNs in the range —100 kms -1 < 
U < 200 km s" 1 of a total bright subsample of 42 PNs; 
however, there is some freedom in the way the angular 
overdensity is removed. In the following, we explore two 
assumptions: (i) using all 12 PNs in the angular over- 
density, but weighting each one by 2/12 = 0.17, and (ii) 
removing the brightest 10 of the 12 PNs in the overden- 
sity. The resulting luminosity functions for both cases 
are plotted in Figure^] they differ only slightly. 

Thus in the following we use the kinematic condition 
— 100 km s -1 < U < 200 km s -1 together with assump- 
tion (ii) above to ensure azimuthal uniformity of the 
bright PNs, as an improved selection criterion for PNs in 
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the main population in NGC 4697. We keep all PNs in 
the intermediate luminosity bin with —100 km s -1 < U < 
200kms _1 , because Fig. [S] showed that their azimuthal 
asymmetry is not large. Also, we have checked that the 
mean angular velocity distribution in this magnitude bin 
after applying the kinematic selection follows approxi- 
mately the sinusoidal variation of the faint PNs, and the 
velocity disperion is approximately constant. The result- 
ing main population sample is identified by their bright- 
ness, distribution and kinematic properties. 

Comparing the luminosity function in Figure ITU1 of this 
main PN population, with the luminosity function of all 
PNs in NGC 4697, we see that the bright cutoff of the 
main population is shifted to fainter magnitudes. We 
note that the bright cutoff could be shifted further to 
fainter magnitudes if some of the kinematically normal 
and azimuthally uniform bright PNs were also part of 
the secondary population, which we cannot tell from the 
present data. 

We can ask now what is the effect, in practice, on 
the PNLF distance determination. The reduced sam- 
ple for the main PN population in NGC 4697 has 214 
objects. After binning these data into 0.2 mag inter- 
vals, we transform the apparent magnitudes m(5007) into 
absolute magn itudes, adopting the extinction correction 
of 0.105 mag (Mendez et al.ll200l|) and assuming differ- 
ent distance moduli, an d we compare the result s with 
the PNLF simulations of lMendez fc Soffnerl lTt9 97L This 
is the same procedure used in lM^ndez"et al . ( 2001) for 
the PNLF distance determination. The comparison is 
shown in Fig. ^2 We conclude that the PNLF distance 
modulus should be increased slightly from 30.1 (the ear- 
lier determination based on the full sample) to 30.2 or 
30.25. For comparison, the x 2 -fit to the same data (blue 
line in Fig. HOf) gives 30.22. This correction would bring 
the PNLF distance modulus in better (but not perfect) 
agreement with the surface brightness fluct uation ( SBF) 
distan ce modulus (30.35 ± 0.2) reported bv lTonrv et al. I 
J200l . 

4. DISCUSSION 

4.1. Can the NGC 4697 PN sample be contaminated? 

Bright PNs play a significant role in the analysis of 
the last section. Is it possible that the brightest PNs 
in NGC 4697 are contaminate d by compact HII re - 
gions such as those obse rved in iGerhard et al] (2002); 
Ryan- Weber et aT| l|2004D ? The recent observations of 
Mendez et al. ( 2005J) have shown that this cannot be so. 
These authors took spectra of 13/42 PNs in our bright 
subsample in NGC 4697 with FORS2@VLT; these have 
no detectable continuum and the line ratios of metal-rich 
PNs. 

The same argument also shows that these bright PNs 
cannot be background emission line galaxies. Moreover, 
Ly alpha emission galaxies would come in at [OIII1 mag- 
nitud es of ni5oo7 <; 26.2 (see Fig. 4 of lAguerri et alJ 
I2005|) . while the bright PNs in NGC 4697 have m 50 o 7 < 
26.2. 

Furthermore, one may wonder whether the bright PNs 
in NGC 4697 might simply be foreground objects which 
could be closer to us and hence brighter than the 'true' 
NGC 4697 PNs on which they would be superposed. 
Given that NGC 4697 is located in the southern ex- 
tension of the Virgo cluster, known to contain an intra- 
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Fig. 10.— Distributions of PN magnitudes in NGC 4697. The 
PNLF for all PNs (black diamonds, das hed black line) is simi- 
lar to the PNLF of lMendez et alJ 120011) (black solid line). The 
red squares and x 2 -fitted dash-dotted line show the PNLF of the 
main population of PNs in the velocity range — lOOkrns -1 < U < 
200 kms -1 , but with the bright PNs in the angular overdensity 
weighted only by 2/12 so that the angular distribution of the bright 
PNs becomes approximately uniform. The blue stars and x 2 -htted 
blue dotted line show the PNLF in the same velocity range, but 
now with the 10 brightest PNs in the angular overdensity removed. 

cluster population of PNs (|Arnaboldi et al J 120021 120041: 
iFeldmeier et al.l2004|) , this possibility deserves to be con- 
sidered. However, the following observational facts show 
that the bright PNs in NGC 4679 are not intracluster 
PNs (ICPNs). (1) NGC 4697 is i n fact closer than the 
Virgo cluster. The PNLF from iMendez et afl l|2001fl 
places it at 10.5 ± IMpc (m-M= 30.1), while the dis- 
tance modulus from the PNLF of M87 is m-M= 30.7±0.1 
ijCiardullo et al J 12002). Even if we discarded the entire 
brightest 0.3 mag of the NGC 4697 PNs, NGC 4697 
would still be at 75% of the distance of M87. (2) The ve- 
locity dispersion of the bright and unrelaxed PNs in NGC 
4697 is ~ 170 km s _1 , but varying azimuthally, while that 
of the fainter main population is 123 km s -1 (Table 
Figs. 0] EJ- Both are much smaller than the veloci ty dis- 
persion of ICPNs in Virgo ijArnaboldi et alJ l2004h The 
radial distribution of the bright population outside the 
incompleteness ellipse is concentrated towards the galaxy 
center and KS compatible (92%) with the radial distri- 
bution of the fainter PNs (see Fig. 0J). Thus also the 
bright PNs in NGC 4697 are bound to the galaxy. (3) 
The surface density of PNs with 1115007 < 26.2 is 0.58 
PNs/arcmin 2 , while the mean surface de nsity of Virgo 
ICPN s is 0.02 PNs/arcmin 2 , much smaller ijAguerri et alJ 
2005). This last argument also rules out significant con- 
tamination by chance superpositions of PNs even closer 
than NGC 4697. 

We conclude that the bright PN population in NGC 
4697 consists of genuine PNs, and that it is dynamically 
bound to the NGC 4697 system. The irregular angular 
distribution and kinematics of these PNs, by comparison 
with the fainter main population, show that they must 
belong to a separate stellar population not yet in dynam- 
ical equilibrium with NGC 4697. 
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Fig. 11— PNLF of NGC 4697 (squares), built from the 
main population PN sample (blue stars in Fig. I1UI of 214 ob- 
jects, binned into 0.2 mag intervals. The apparent magnitudes 
m(5007) have been transformed into absolute magnitudes M (5007) 
by adopting an extinction correction of 0.105 mag and distance 
moduli indicated in each plot. The three lines are PNLF simula- 
tions (Mcndcz & Soffncr 1997). For a distance modulus 30.15 the 
brightest PNs are a bit too weak, therefore the distance must be 
increased. But 30.35 is clearly excessive. The best fit is for 30.2 or 
30.25. 

4.2. Comparison with other galaxies 

PN samples as large as the one for NGC 4697 are still 
the exception. Yet to undertake the analysis described in 
Section |3 it was crucial to work wit h a complete sample 
of some 300 PNs. The PN.S project Ipouelas et al.l2002ft 
may lead to complete samples of simil ar size but as of now 
the t ypical sample sizes are < 200 (|Romanowskv et al.l 
2003) and their magnitude distributions and complete- 
ness have not been studied. The only other elliptical 
gal axy with a comparable (even larger) sample is Cen 
A l|Hui et al.lll995t IPeng. Ford fe Freemanll2004j) . A de- 
tailed analysis of the distribution of Cen A PNs in the 
magni tude-velocity plane is still pending, but iHui et aTl 
( 1993) analysed the PNLF in Cen A as a function of ra- 
dius, based on narrow-band surveys. They concluded 
that no population effect on the PNLF bright cut-off 
could be seen, suggesting that the 0.3 mag difference 
between their main and outer halo samples was due to 
filter transmi ssion uncertainties. Analysis of the large 
sample from IPeng. Ford fc Freemanl l)2004|) will clarify 
whether the PNs in Cen A are consistent with one or 
more subpopulations. This will be particularly inter- 
esting because Cen A is believed to be the remnant of 
a galaxy merger, so one might expect PNs from both 



the older stellar populations of the progenitors as well 
as from the stars formed in the subsequent interaction 
between them. 

Without large kinematic samples, searches for PN sub- 
populations must be based on the PN luminosity dis- 
tributions. The work of Jacoby, Ciardullo and collab- 
orators cited in the Introduction has shown that the 
PNLF is remarkably uniform. However, there are ex- 
ceptions: we recall that in the halo of M84 there exists 
a small population of overluminous PNs whose cutoff is 
0.3 mag brighter than that of the main M84 population, 
but which appear nonetheless bound to the halo of M84 
ijArnaboldi et alJl200l lAguerri et alJl2005j) . They must 
therefore be intrinsically bright, due to some stellar pop- 
ulation difference. In M87, the only overluminous PNs 
projected onto the galaxy for whi ch velocities have been 
measured ijArnaboldi et alJ l2004f) . have very large rela- 
tive velocities with respect to M87. This is most natu- 
rally explained if these PNs have fallen into the deep po- 
tential well of M87 from far out in the cluster; this would 
again imply an intrinsic population difference. However, 
a larger kinematic sample in M87 is required to be cer- 
tain. 

The existence of the bright PN subpopulation in NGC 
4697 implies some uncertainty in the PNLF cut-off lu- 
minosity for this galaxy. Depending on whether or not 
the azimuthally symmetric brightest PNs are part of 
the main population, the cut-off luminosity of the main 
PN population in NGC 4697 is fainter than that of the 
whole population by > 0.15 mag (Figs. lllfl . While 
this is consistent with the M84 result, is it also con- 
sistent with the system atic studies of PNLF distances? 
ICiardullo et alJ (|2002f) have compared the PNLF and 
SBF distances, finding a distribution of residuals with 
a systematic offset by ~ 0.3 mag, which they suggested 
is due to internal extinction effects, and with a FWHM of 
~ 0.5 mag. The offset has in the mean-time been reduced 
to ~ 0.1 5 mag, follow i ng a r evision of the SBF distance 
scale by Uensen et alJ l)2003|) . The width of this distri- 
bution is consistent with their determination of the ob- 
servational errors in both methods. However, the offset 
we have determined for NGC 4 697 is also consistent w ith 
the distribution of residuals in ICiardullo et ail (J2002). 

4.3. Origin of the PN population difference: a 
secondary, younger stellar population? 

We have shown that a large fraction of the bright 
PNs in NGC 4697 belong to a secondary, dynamically 
young stellar population that is not well-mixed in the 
gravitation al potential of the galaxy. Late infall of tidal 
structures (jZezas et alJ 2003) or a merger with a smaller 
galaxy some time ago would be natural ways to add such 
an unmixed stellar component to NGC 4697. What phys- 
ical population difference is correlated with this dynam- 
ic al youth? 

Mendcz et al. 



(|2005|) show from their spectroscopic 
data for 13 bright PNs that these have near-solar metal- 
licities. Of these 13 bright PNs, 6 are inside the incom- 
pleteness ellipse, one has no measu red velocity, and 6 be- 
long to our secondary population. iMendez' eta!] (2005) 
also use long-slit spectroscopy to show that the metallic- 
ity of the integrated stellar population within one effec- 
tive radius has solar or higher metallicities. These ob- 
servations make it unlikely that metallicity is the main 
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factor responsible for the different magnitude distribu- 
tions of the main and secondary PN populations in NGC 
4697. 

Thus the more likely driver would appear to b e an 
age difference, as suggested bv lMarigo et all l)2004D and 
as might generally be expected in an accretion event. 
Based on their result that the distribution of X-ray point 
sources in NGC 4 697 does not follow the stellar light, 
IZezas et"aLl (I2003D have argued that this is because these 
sources were formed several 10 8 years ago in tidal tails 
that are now falling back onto the galaxy. Note, however, 
that the integrated light in NGC 4697 sho ws no evidence 
of yo ung stars with mean age < 7 Gyr ijMendez et alJ 
2005), so this younger component could not be lumi- 
nous enough to contaminate the integrated light to the 
level measured. Also note that the observed increase 
of extinction in the PN envelope with PN core mass 
more than compensates for the increase of core luminos- 
ity with core mass, for brig ht PNe in local group galaxies 
ijCiardullo fc Jacobvll999|) . so that stars with ages below 
1 Gyr may not reach the [OIII] luminosity at the PNLF 
cutoff. A secondary stellar population younger than 1 
Gyr is ther efore unlikely as well. 

Recently, Ciard ullo et alJ l)2005f) have argued that the 
brightest PNs in the PNLF must have core masses of 
> 0.6 Mq, corresponding to main sequence masses of 
~ 2.2 M . They argue further that for such high-mass 
objects to occur in elliptical galaxies, these early-type 
galaxies would either have to contain a small, smoothly 
distributed component of young (< 1 Gyr age) stars, or 
more likely, that the bright PNs in these systems have 
evolved from blue straggler stars created through binary 
evolution. Their blue straggler model, due to the as- 
sumption of a fixed distribution of primary-to-secondary 
mass ratios for the initial binaries, predicts that older 
stellar populations produce fewer bright PNs per unit 
luminosity, as is observed, because the number of bi- 
nary stars in a stellar population that can coalesce to 
~ 2.2 Mq blue stragglers decreases with time. 

If correct, this blue straggler model could also ex- 
plain how the secondary population we found in NGC 
4697 can contain a large fraction of the brightest PNs in 
this galaxy, provided that the stellar population corre- 
sponding to this secondary PN population is appreciably 
younger than the main stellar popu lation, whose age i s 
~ 9 Gyr from optical spectroscopy ijTrager et alJ l20001. 
At the same time, this secondary stellar population must 
not be so young to violate either the constraints from the 
optical colours or from the envelope absorption - PN core 
mass correlation, i.e., must be older than ~ 1 Gyr. We 
can give an estimate for the effect of such an intermedi- 
ate age population on the optical colours as follows. The 
secondary subpopulation traced by the bright and pre- 
dominantly counterrotating PNs may contain ~ 20% of 
all PNs in the Complete sample for NGC 4697. A stellar 
population as blue as the bulge of M31 has a luminosity- 
specific PN density per unit Lb up to 5 times higher than 
the populations characteristic for old elliptical galaxies 
l|Hui et all 11991 iCiardullo et al.ll2005|) . Thus the sub- 
population corresponding to the secondary PN popula- 
tion in NGC 4697 would be expected to contain > 5% 
of the blue luminosity of NGC 4697, spread over a large 
fraction of at least the E image. To detect this we need 
deep and accurate photometry. 



The unmixed spatial and velocity distributions of the 
secondary PN population in NGC 4697 show that that 
this population is dynamically young, i.e, has not had 
time to phase-mix and come to dynamical equilibrium 
in the gravitational potential of NGC 4697. It may well 
be associated with tidal structures that were formed in 
a merger or accretion event ^1 — 2 Gyr ago, and that 
have now fallen back onto the galaxy, or be associated 
with a more recent merger /accret i on wit h a red galaxy 
such as described in Ivan Dokkum I ((2005) . In a Universe 
in which structures form hierarchically, such secondary 
stellar populations might be quite common in ellipitical 
galaxies, but they would be difficult to see. The present 
work shows that studying their PN populations is one 
promising approach of looking for such secondary popu- 
lations. However, large PN samples are required; most 
existing PN studies of early-type galaxies do not have the 
statistics for such an investigation. Moreover, in only a 
fraction of cases may there be enough asymmetry signal 
to detect with a few hundred PNs. 

5. CONCLUSIONS AND IMPLICATIONS 

We have analysed the magnitudes, kinematics and po- 
sitions of a complete sa mple of 320 PNs in t he elliptical 
galaxy NGC 4697 from IMendez et all pOOl^ . This data 
set is large enough for drawing statistically significant 
conclusions, and it does not suffer from detection incom- 
pletenesses in either magnitudes or radial velocities. We 
know of no systematic effects in the data that could ex- 
plain our results. Our main conclusions are: 

1. Bright and faint PNs in NGC 4697 have signifi- 
cantly different radial velocity distributions. The 
mean velocities of the faint and bright subsamples 
(co-rotating and ~ 0, respectively) and their veloc- 
ity dispersions are different, with 94% and 99.3% 
confidence. Thus the PNs in NGC 4697 do not con- 
stitute a single population that is a fair tracer of 
the distribution of all stars. 

2. The luminosity functions of the extreme counter- 
rotating subsample (U < — 100 km s -1 ) and of the 
main population (defined by — 100 km s -1 < U < 
200kms _1 ) are different with 99.7% confidence. 
The PNLF is therefore not universal. 

3. Based on this, we suggest that there exist (at least) 
two PN populations in this galaxy. The secondary 
PN population in NGC 4697 is prominent in a sub- 
sample of counter-rotating PNs brighter than 26.2. 
The luminosity function of the entire extreme coun- 
terrotating sample may be a first approximation to 
the luminosity function of the secondary popula- 
tion. 

4. The spatial distribution of bright PNs with 
to(5007) < 26.2 is different from that of the faint 
PNs. The bright PNs do not follow the azimuthal 
distribution of the optical light, show a left-right 
asymmetry, and have a positive mean radial veloc- 
ity on both sides of the galaxy major axis, but zero 
velocity and larger dispersion on the minor axis. 
They are not in dynamical equilibrium in the po- 
tential of the galaxy. The fainter population has 
rotation properties more similar to the absorption 
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line velocities, with azimuthally constant disper- 
sion. 

5. Using both their kinematics and angular distribu- 
tion, we can estimate a lower limit to the statistical 
fraction of bright PNs in the secondary population. 
Based on this we estimate that the bright cutoff of 
the main population is uncertain by ~ 0.15 mag. 

Our results have two main implications for the use of 
PNs in extragalactic astronomy. First, for distance de- 
terminations with the PNLF, it may be important to 
understand how uniform the PN populations in the tar- 
get galaxies are. From our analysis in NGC 4697 we esti- 
mate that unrecognized subpopulations of PNs in smaller 
samples than that in NGC 4697 may lead to variations of 
~ 0.15 mag in the bright cutoff. This would correspond 
to distance errors of some 10%, which, although a minor 
effect, could be significant in some cases. It will be neces- 
sary to verify how frequently such subpopulations occur 
in elliptical galaxies. We recall that also in the halo of 
M84 there exists a small population of overluminous PNs 
whose cutoff is 0.3 mag brighter than that of the main 
M84 population, b ut which appear nonetheless bound 
to th e halo of M84 ijArnaboldi et al.ll2004l: lAeuerri et all 
2005). 

The second implication concerns the use of PNs as trac- 
ers for the angular momentum and gravitational poten- 



tials of elliptical galaxies. Our analysis has shown that 
in NGC 4697 the bright PNs do not trace the distribu- 
tion and kinematics of stars and are not in dynamical 
equilibrium in the gravitational potential of the galaxy. 
The fainter PNs look more regularly distributed but may 
also contain a fraction of PNs that belongs to this out- 
of-equilibrium population. Clearly therefore, mass deter- 
minations based on PN kinematics will in future require 
careful study of the PN samples being used, not only to 
verify that these PNs are in dynamical equilibrium, but 
also to test for different dynamical components. Even 
if in equilibrium, a younger population of stars may be 
more flattened or have a steeper fall-off than the main 
body of the elliptical galaxy. If the PNs from this popu- 
lation are indeed somewhat brighter than the main pop- 
ulation, one can recognize such differences from their 
lower velocity dispersion or different radial density pro- 
file. However, deep observations and large PN samples 
will be required. 
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